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Abstract 

This paper includes results centered around three topics, all of them related with the nonlinear 
stability of equilibria in Poisson dynamical systems. Firstly, we prove an energy- Casimir type 
sufficient condition for stability that uses functions that are not necessarily conserved by the flow 
and that takes into account certain asymptotically stable behavior that may occur in the Poisson 
category. This method is adapted to Poisson systems obtained via a reduction procedure and we 
show in examples that the kind of stability that we propose is appropriate when dealing with the 
stability of the equilibria of some constrained systems. Finally, we discuss two situations in which 
the use of continuous Casimir functions in stability studies is equivalent to the topological stability 
methods introduced by Patrick et al. |Paal02| . 

1 Introduction 

The use of the conserved quantities of a Hamiltonian flow in the study of the stability of its solutions is 
a venerable topic that goes back to Lagrange and Dirichlet. In the past decades these ideas have been 
adapted to various setups: equilibria in Poisson systems A66, Hoal85, Paal02 , relative equilibria Pa92 
LS98, Or98 , OrRa99 , Paal02 and periodic and relative periodic orbits OrRa99a, OrRa99b of symmetric 
Hamiltonian systems, relative equilibria of symmetric Lagrangian systems SLM91 , and symmetric 
constrained systems Zcal98 , to list a few. All these results provide sufficient conditions for the solution 
in question to be stable. 

In this paper we will focus on the stability of the equilibria of Poisson dynamical systems. The 
main tools that one finds in the literature concerning this case are the energy-Casimir method and the 
topological stability methods introduced in Paal02 . The energy-Casimir method consists in finding 
a combination of conserved quantities by the Hamiltonian flow, typically the energy and the Casimir 
functions, that exhibits a critical point at the equilibrium with definite Hessian. Since the dynamics of 
the system is confined to the level sets of this combination and, by the Morse Lemma, in a coordinate 
chart about the equilibrium these level sets are diffeomorphic to spheres centered at the equilibrium, 
stability follows. The topological methods in Paal02 rely on a much more subtle confinement of the 
flow that takes advantage not only of its conservation laws but also of the topological properties of the 
foliation of the Poisson manifold by its symplcctic leaves. 

Energy confinement is a very important tool in the symplectic Hamiltonian context due to the 
absence of asymptotically stable behavior. Energy methods are, to this day, the only general way to 
prove stability in more than two degrees of freedom. The conservation of the phase space volume by 
the flow imposed by Liouville's theorem does not necessarily hold in the Poisson category. The first 
main result of this paper, contained in Theorem l2.5l adapts the standard energy-Casimir method to this 
situation to allow the use of functions that are not necessarily conserved by the flow but that can still 
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be used to conclude a certain kind of asymptotic stability via the standard Lyapunov stability theorem. 
This newly introduced notion of stability implies the standard Lyapunov stability and will be referred 
to as weak asymptotic stability. Theorem 12.51 improves a previous version of the energy-Casimir 
method (see |Or98] or Corollary 4.11 in OrRa99b ) where the conserved quantities confining the flow 
are also used to shrink the space on which one checks the definiteness of the Hessian. Theorem 12.51 
shows that any conserved quantity can be used to shrink this space even when that conserved quantity 
is not involved in the confinement of the flow. 

Theorem l2.1()l is the second main result of this paper. It adapts the stability condition in Theorem l2.5l 
to equilibria of Poisson systems obtained by a certain reduction procedure that uses ideals in the Poisson 
algebra of the functions on the manifold. Our interest is twofold. First, there are some constrained 
mechanical systems that can be described by reducing in this sense a bigger (unconstrained) system. 
Moreover, the weakened kind of stability that Theorem 12.101 allows us to conclude, coincides with the 
physically relevant notion of stability in those situations, that is, the one that describes the system when 
subjected to perturbations compatible with the constraints. We illustrate this point with a couple of 
examples in Section|21 a light Chaplygin sleigh on a cylinder and two coupled spinning wheels. Second, 
there are cases when there are not enough conserved quantities to apply Theorem 12.51 but . nevertheless, 
the system can be reduced around the equilibrium and then the reduced system has enough conserved 
quantities to use the theorem. Theorem 12.101 explains the meaning of having this reduced kind of 
stability. In particular, it shows the role of sub-Casimir functions in stability computations. 

The last section of the paper is dedicated to the study of the relation between the topological stability 
methods in |Paal02| with a generalized version of the energy-Casimir method that we propose in the 
text based on the use of local continuous Casimir functions of the Poisson manifold. To be more explicit, 
the stability criteria in |Paal02| are stated in terms of a set that, roughly speaking, measures how far 
the space of symplcctic leaves of a Poisson manifold is from being a Hausdorff topological space. The 
general question that we try to answer is under what circumstances this set can be characterized as the 
intersection of level sets of local continuous Casimirs. Since this is not true in general, we provide two 
sufficient conditions that are related to certain idempotency of the set in Paal02 and to the possibility 
of separating regular symplectic leaves by using continuous Casimirs. The natural category where 
these questions are posed is that of generalized foliated manifolds; this is the context in which we have 
formulated the main results in this section and where we have obtained the Poisson case as a byproduct, 
considering it as a manifold foliated by its symplectic leaves. 

Before we start with the core of the paper we quickly review in a few paragraphs the basic notions 
and terminology of Poisson manifolds and generalized foliations that we will use throughout the paper. 
The expert can safely skip the rest of this section. 

Poisson systems. Let P be a smooth manifold and let C°°(P) be the algebra of smooth functions on 
P. A Poisson structure on P is a bilinear map {•, •} : C°°(P) x C°°(P) — ► C°°(P) that defines a Lie 
algebra structure on C°°(P) and that is a derivation on each entry. The derivation property allows us 
to assign to each function F £ C°°(P) a vector field Xp £ X(P) via the equality 

X H [F] := {F, H} for every F £ C°°(P). 

The vector field Xh £ X(P) is called the Hamiltonian vector field associated to the Hamiltonian 
function H. The derivation property of the Poisson bracket also implies that for any two functions 
F, G £ C°°(P), the value of the bracket {F, G}(z) at an arbitrary point z £ P depends on F only 
through dF(z) which allows us to define a contravariant antisymmetric two-tensor B £ A 2 (P) by 

B(z)(a z , (3 Z ) = {F, G}(z), 

where dF(z) — a z £ T*P and dG(z) — (3 Z £ T*P. This tensor is called the Poisson tensor of M, The 
vector bundle map 

B » . T * p Tp naturally assoc i a ted to B is defined by B(z)(a z , f3 z ) = (a z , B^(p z )}. 
Its range £ := B^(T*P) C TP is called the characteristic distribution of the Poisson manifold 
(P, {-,-}). Its fiber at z £ P is hence given by S z = {X H (z) | H £ C°°(P)}. Th e distr i bution £ is 
a smooth generalized distribution which is always integrable in the sense of Stefan |St74al ISt74b| and 
Sussmann |Su73| . Its maximal integral submanifolds {£} are symplectic and are called the symplectic 



Ortega, Planas-Bielsa, and Ratiu: Asymptotic and Lyapunov stability of Poisson equilibria 



3 



leaves of (P, {•, •}). The symplectic form uj£ on the leaf C is uniquely characterized by the identity 

lu c (z) (X F {z),X G (z)) := {F, G}(z) for any F, G £ C°°{P) and for any z e £.. 

Since the symplectic leaves of (P, {•,•}) are the maximal integral leaves of a generalized distribution, 
they form a generalized foliation in the sense of |Daz85| . This implies the the existence of a chart 
([/, ip : U — > M. m ) around any point z € P such that if £ 2 is the symplectic leaf containing z then there 
is a countable subset A C M m_ ", with m = dimP and n = dim£ 2 , such that 

<p{U HC z ) = {ye tp(U) I (y" +1 , . . . , y m ) £ A}. (1.1) 

Such a chart (U, (p) is called a foliation chart for the generalized symplectic foliation of P around 
the point z. A connected component of U D C z is called a plaque of the foliation chart (U,ip). The 
point z is said to be regular if the neighborhood U can be shrunk so that all the leaves that it 
intersects have all the same dimension. In that case, the plaques coincide with the points of the form 
(y 1 , . . . , y n , 2/o +1 j ■ ■ ■ , y™) £ <f(U) with (yQ +1 , . . . , y™) constant. A leaf consisting of regular points is 
said to be regular and singular otherwise. The set of regular points of a generalized smooth foliation 
is open and dense. 

Some of the results proved in this paper will be first given in the category of foliated manifolds. The 
corresponding results in the context of Poisson manifolds are then obtained as corollaries. 

Casimirs, local Casimirs, and first integrals of foliations. A function on a foliated manifold that 
is constant on the leaves is called a first integral of the foliation. When we consider the particular case 
of a Poisson manifold, the elements in the center of the Poisson algebra (C°°(P), {•,•}), also called the 
Casimir functions, are first integrals of the foliation of P by its symplectic leaves. A local Casimir 
at the point z £ P is a function C £ C°°(U Z ) for some open neighborhood U z C P of z such that it is a 
Casimir of the Poisson manifold (U z , {•, -}u z ) where the bracket {•, -}u r is the restriction of the bracket 
{•, •} on P to U z . 

In general, non-trivial global Casimir functions may not exist. On the other hand, local Casimirs 
are always available in the neighborhood of a regular point. Indeed, if we think of the Poisson manifold 
(P, {•,•}) as a foliated space by its symplectic leaves, the expression 1)1.1(1 allows us to find a chart 
(U, ip : U — > W n ) around the regular point where the plaques of the symplectic foliation are the points 
of the form (y 1 , . . . , y n , j/q +1 , . . . , y™) £ <p{U) with (j/q +1 , . . . , y™) constant. The functions that depend 
on the last m — n coordinates are local Casimir functions of (P, {•, •}) around z. 

Quasi-Poisson submanifolds and sub-Casimirs. An embedded submanifold S of P which is Poisson 
in its own right and is such that the inclusion i : S P is canonical is called a Poisson submanifold 
of P. The Poisson structure on S is uniquely determined by the condition that the inclusion be canonical, 
that is, there is no other Poisson structure on S relative to which the inclusion is canonical. 

It turns out that in this paper we need a slightly weaker condition. An embedded submanifold S 
of P (without any Poisson structure on it) such that P"(s) (T*P) C T S S for any s £ S is called a 
quasi-Poisson submanifold of P. Every Poisson submanifold is quasi-Poisson but the converse is 
not true. As a corollary to the main theorem in MaRa86 , one can easily conclude that if S is a quasi- 
Poisson submanifold of P, then there is a unique Poisson structure {•, -}s on S with respect to which 
the inclusion S P is a Poisson map, that is, there is a unique induced Poisson structure on S making 
it into a Poisson submanifold of P. The Poisson bracket {•, -}s is defined by {/, g}s(s) '■= {F,G}(s) 
where P, G £ C°°(P) are arbitrary local extensions of /, g £ C°°(S) around the point s £ 5; this means 
that there is an open neighborhood U of s in P such that f\snu = F\snu and g\snu — G\snu- 

Thus, it is possible that the quasi-Poisson submanifold S of P has its own Poisson structure (that 
is given a priori) but it is not the one induced by the Poisson structure of P. For a discussion of these 
issues see |OrRa03| . sections 4.1.21 - 4.1.23. 

Let c £ C°°(S) be a Casimir function for the Poisson manifold (S, {-, -}s)- Any extension C £ C°°(P) 
of c will be called a sub- Casimir of (C°°(P), {•,•})• 

Here is an example of the construction just described. Take some Casimir functions C\, . . . ,Ck £ 
C°°(P) of (P, {•,•}) and assume that a certain common level set S of these Casimirs is an embedded 
submanifold of P. It is easy to check that B*(s) (T*P) C T S S for any s £ S and hence S carries a unique 
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Poisson bracket {•, -}s) such that (S, {■, -}s) is a Poisson manifold with its own Casimir functions that 
extend to sub-Casimirs on P. 

2 Stability in Poisson systems 

In this section we use some aspects of the geometry of Poisson manifolds to study the stability of the 
equilibria of Hamiltonian vector fields. 

Let M be a manifold, X £ X(M ) a vector field, F t the flow of X, and m e £ M an equilibrium of X, 
that is, X(m e ) — or, equivalently, F t (m e ) — m e for all t £ R. Recall that m e is stable, or Lyapunov 
stable, if for any open neighborhood J7 of m e in M there is an open neighborhood V C U of m e such 
that F t (m) £ J7 for any m £ V and for any £ > 0. The equilibrium m e is asymptotically stable if 
there is a neighborhood V 7 of m e such that -Ft(V) C F S (V) whenever t > s and lim F t (V) = m e , that 

i — >oo 

is, for any neighborhood W of m e there is a T > such that F t (V) C W ii t > T. If only the first 
condition holds, that is, F t (V) C F S (V) whenever t > s, we say that m e is weakly asymptotically 
stable. Note that 

asymptotic stability => weak asymptotic stability => Lyapunov stability. 

Asymptotic stability cannot occur in symplectic Hamiltonian systems due to Liouville's theorem ; only 
Lyapunov stability is allowed. In the Poisson category, equilibria lying in trivial symplectic leaves 
may be asymptotically stable. However, if the symplectic leaf that contains the equilibrium is at least 
two-dimensional, weak asymptotic stability is the most we can hope for. 

The linearization of X at the equilibrium point m e is the linear map L : T mis M — > T m ^M defined 
by L(v) := Jj| t _ (T me F t (v)) where F t is the flow of X and v £ T m<i M is arbitrary. As is well known, 
the study of the spectrum of the linear map L gives relevant information about the stability of the 
equilibrium m e . The equilibrium m e £ M is linearly stable (respectively unstable) if the origin is a 
stable (respectively unstable) equilibrium for the linear dynamical system on T mr M defined by L. The 
equilibrium m e is spectrally stable (respectively unstable) if the spectrum of the linear map L lies 
in the (strict) left-half plane or on the imaginary axis (respectively at least one eigenvalue has strictly 
positive real part). Lyapunov and linear stability imply spectral stability. If all the eigenvalues of L have 
strictly negative real part, that is, they lie in the (strict) left-half plane, the system is asymptotically 
stable. 

2.1 Linearization of Poisson dynamical systems and linear stability 

Consider a Hamiltonian vector field Xh on the Poisson manifold (P, {•,•}), let z e £ P be an equilibrium 
of Xh, and L : T Zr P — > T Zc P the linearization of Xh at z e . If z e is regular (in particular, when P is a 
symplectic manifold) there are restrictions on the eigenvalues of L that do not allow us to conclude the 
Lyapunov stability of z e from its spectral stability (see, for instance, Theorem 3.1.17 in |AM78| 1. As 
will be shown below, this restriction disappears, in general, for equilibria lying on singular symplectic 
leaves. 

In order to present the following lemma, whose proof is a straightforward computation, we recall that 
there exists a chart (U, tp) around any point z £ P in the 2n + r dimensional Poisson manifold (P, {•, •}) 
such that (p(z) = and that the associated local coordinates, denoted by (q , . . . , q n ,pi, . . . ,p n , z , . . . , z r ), 
satisfy {q l , qi} = {pi,Pj} = {q l , z k } = {pi, z k } = and {q l ,Pj} = 5), for all i, j, k such that 1 < i,j < n, 
1 < k < r. For all such that k,l, 1 < k,l < r, the Poisson bracket {z k ,z 1 } is a function of the local 
coordinates z , . . . , z r exclusively and vanishes at z. Hence, the restriction of the bracket {•, •} to the 
coordinates z 1 , . . . , z r induces a Poisson structure on an open neighborhood V of the origin in W whose 
Poisson tensor will be denoted by 1Z £ A 2 (V). This Poisson structure on V is called the transverse 
Poisson structure of (P, {•, •}) at z and is unique up to Poisson isomorphisms. The coordinates of 
the local chart that we just described are called Darboux- Weinstein coordinates We83 . 
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Lemma 2.1 Let z e be an equilibrium of the Hamiltonian dynamical system on the Poisson manifold 
(P, {•,•}) and let (q, p, z) be a Darboux-Weinstein chart around z. Denote by x := (q, p) and by J the 
n x n square matrix given by 

HI o)- 

The linearization L of Xh at the equilibrium z e in the coordinates (x, z) takes the form 
where 

p— 1 p=l p=l 

Proof. The result is obtained by differentiating the expression of the Hamiltonian vector field at the 
equilibrium in Darboux-Weinstein coordinates and by taking into account that the matrix J is constant, 
that 72.(0) is zero, and that 1Z depends only on the z variables. ■ 

We now use (|2.1|l to give a characterization of the structure of the eigenvalues of the linearized vector 
field L in the Poisson context. The proof of the following proposition is a straightforward computation. 

Proposition 2.2 In the situation described in the previous lemma denote by {Ai, . . . , \2n} the eigenval- 
ues of the infinitesimally symplectic matrix S, counted with their multiplicities, and let {ui, ■ ■ ■ , i«2n} be 
a basis of corresponding eigenvectors. Assume that the matrix V is diagonalizable, let {fM\, . . . , /i r } be its 
eigenvalues counted with their multiplicities, and {y\, . . . ,v r } a basis of eigenvectors. Then the matrix 
L has eigenvalues {Ai, . . . , A2„, fii, ■ ■ . , A* r }. If for any eigenvalue fij we have that (S — [ijl)^ 1 Qvj is 
not empty then L is diagonalizable with corresponding basis of eigenvectors 

{(ui,0), . . . (u 2n ,0), (-wi,vi), (-w r ,v r )}, 

where Wj G (<S — HjI)~ 1 Qvj , j = 1, ■ ■ • ,T are arbitrary but subjected to the condition that if Vj = Vk 
then (ujj,Vj) and (wk,Vk) are chosen to be linearly independent. 

The eigenvalues {Ai, . . . , A2 n } satisfy the symplectic eigenvalue theorem since S is infinitesimally 
symplectic. However, the eigenvalues {fJ>i, ■ ■ ■ , /i r } may lie, in principle, anywhere in the complex plane. 
Hence Poisson dynamical systems may exhibit asymptotic behavior. There are three specific situations 
that should be singled out: 

• None of the eigenvalues of V coincides with one of the eigenvalue of S. In this case the matrices 
[S — fj,jl), 1 < j < r, are invertible and the whole linear system L is diagonalizable. 

• Hi = Xj for some i, j but (S — HiI)~ 1 Qvi is not empty. Then there is a passing of eigenvalues but 
they do not interact in the sense that they correspond to different blocks in the linearized system. 
We will call this situation uncoupled passing. 

• If in the previous case (<S — ^iI)~ 1 Qvi is empty then the linear system is not diagonalizable 
anymore and the passing of eigenvalues mixes blocks of the infinitesimally symplectic part and the 
transversal one. We will call this situation coupled passing. 

With these remarks in mind, we get the following. 

Proposition 2.3 Let (P, {•,•}, H) be a Poisson dynamical system and z e £ P an equilibrium point of 
Xh . If the linearization L of Xh at z e exhibits a coupled passing then the system is linearly unstable. 

Proof The existence of a coupled passing implies the occurrence in L of a nondiagonal block in its 
Jordan canonical form. The flow of the linear dynamical system induced by L, when restricted to the 
space generated by the associated Jordan basis, exhibits an unstable behavior and the result follows. ■ 
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Corollary 2.4 Consider the linearization L of a Poisson dynamical system (P, {•,•}, H) around an 
equilibrium z e £ P lying on a regular symplectic leaf C. Let {X%, . . . , X^n} be the eigenvalues of the 
infinitesimally symplectic block S. Then 

(i) v = o. 

(ii) The vectors u € T Za P that satisfy Lu — Xu for some X ^ lie in T z<! £. In particular, the unstable 
directions of L are tangent to the symplectic leaf of P that contains the equilibrium. 

(iii) If S~ 1 Qvj is not empty for any Vj as in Provosition \2.2\ then is the only eigenvalue in addition 
to {Ai, . . . , X 2n }. 

Proof. The first part follows from the expression for V provided in Lemma I2"T1 and from the fact that 
TZ = in an open neighborhood of z e that contains only regular points. The unstable directions are the 
vectors in the eigenspaces corresponding to strictly positive eigenvalues. Then the points (ii) and (iii) 
follow from the expression of L in Lemma \2 . II using that on the set of regular points 71 = V = 0. ■ 

2.2 Non linear stability in Poisson dynamical systems 

As noted in the previous subsection, the array of linear tools available to conclude nonlinear stability 
of equilibria of a Poisson dynamical system is very limited. In this section we will state a result that 
provides a sufficient condition for such equilibria to be Lyapunov or weakly asymptotically stable. This 
result is based on the use of local Casimirs and conserved quantities of the dynamical system in question 
and is related to the classical energetics methods (also called Dirichlet criteria) in |A66I IPaal02| . Our 
approach builds on an improvement of the classical result in |A66| that was carried out in [Or98| (see 
Corollary 4.11 in jOrRa99bp . 

The proof of our main result will be based on a classical result of Lyapunov that states that if 
m e G M is an equilibrium of the vector field X G X(M ) with flow F t and there exists a positive function 
L G C°°(U) around m e , with U an open neighborhood of m e , such that L(m) := | . L(F t (m)) < 0, 
for any m G U\ {m e }, then m e is a Lyapunov stable equilibrium. We recall that a function / G C°°(M) 
is said to be positive around m e G AI if f(m e ) = and there is an open neighborhood U me of m e such 
that f(m) > 0, for all m G U m ,, \ {m e }. If L(m) < for all m G U mc \ {m e }, then m e is asymptotically 
stable. See e.g. Theorem 1, Chapter 9, §3 in HS74 for a proof of these statements; the infinite 
dimensional versions of these assertions can be found in Theorems 4.3.11 and 4.3.12 of |AMR88j . Any 
positive function L in the statement of Lyapunov's theorem is usually called a Lyapunov function. 
Its construction for specific dynamical systems is by itself a very active research subject. 

In the case of Hamiltonian mechanics, the Hamiltonian and the Casimirs of the Poisson phase space 
are natural candidates to be used in Lyapunov's theorem. If, additionally, the system has a symmetry to 
which one can associate a momentum map, its components are conserved quantities that sometimes can 
be used for the same purpose. The use of all conserved quantities of a dynamical system in the study 
of the stability of equilibria to form Lyapunov functions is known under the name of energy-momentum 
methods. However, it should be noted that, apart from conserved quantities, Lyapunov's theorem can 
be applied with the more general class of functions whose time derivative is strictly negative. The 
existence of these functions implies the asymptotic stability of the equilibrium in question. In the 
symplectic context this is impossible. This behavior, allowed only for Poisson Hamiltonian systems, is 
used in the main theorem of this subsection. 

In the sequel we will use the following notation. Let P be a smooth manifold, / G C°°(P) a smooth 
function, z e G P a critical point of / (that is, d/(z e ) = 0), and U an open neighborhood of z e . The 
Hessian of / at the critical point z e is the symmetric bilinear form d 2 f(z e ) : T Ze P x T Ze P — > M given 
by d 2 f(z e )(v, w) := where v,w G T Zc P and W G X(U) is an arbitrary extension of w to a 

vector field on U. The fact that z e is a critical point of / ensures that this definition is independent of 
the extension W of w. 

Theorem 2.5 Let (P, {•, •}, H) be a Poisson dynamical system, z e an equilibrium point of the Hamilto- 
nian vector field Xh, and C\, . . . ,C'k ' P — ■> R conserved quantities, that is, {Ci, H} = 0, i € {1, . . . , k}. 
Let F : P —>M. be a function satisfying at least one of the following two conditions 
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(i) {{F-F{z e )) 2 , H}(z) < for all z G P\{z e } and {F-F(z e ), H}(y) < for all the points y G P\{z e } 

such that {{F - F(z e )) 2 , H}(y) = 0. 

(ii) {(F — F(z e )) 2 , H}(z) < for all z G P \ {z e } and F ~ F(z e ) is positive on P around z e . 
Assume that there exist constants {Xq, Ai, ■ ■ ■ , fi} such that 

d(X H + AiCi + . . . + X k C k + nF)(ze) = 

and the quadratic form 

d 2 (X H + X 1 C 1 + ... + X k C k + f iF)\ WxW {z e ) (2.2) 

is positive definite, where 

W = kerdH(z e ) n kerdCi(^ e ) n . . . n kerdC fe (z e ). 

Then z e is a weakly asymptotically stable equilibrium (and hence Lyapunov stable). If the inequality 
{(F - F{z e )) 2 ,H}(z) < is strict for every z G P\{z e } then z e is asymptotically stable (this can only 
happen if the symplectic leaf that contains the equilibrium is trivial). 

Proof Consider the functions G C°°(P) denned by 

k 

h(z) := J2 (^Cj(z) + »F(z)) - (Aj-QCO + fiF(z e )) , 

Hz) -E^ ((^w - c iW) J + ™ - F (^)) 2 ) ■ 

3=0 

where we denote Cq := H . Notice that h(z e ) — and that, by hypothesis, dh(z e ) — which implies that 
d 2 /i(z e ) is well defined. Moreover, hypothesis (|2.2|) is equivalent to the statement that d 2 li(z e )\wxw 
is positive definite. Additionally, l2(z&) = 0, dl2(z e ) = 0, and hence d 2 l2(z e ) is well defined. A 
straightforward computation shows that d 2 l 2 {z e ) is positive semidefinite with kernel equal to the space 
W. A result due to Patrick (see |Pa92j ) shows that in these circumstances there exists a constant r > 
such that for any e G (0,r] the Hessian d 2 (eli + l2)(z e ) is positive definite. 

Let L e := el\ + 1%. The positive definiteness of d 2 L e (z e ) implies that L e is a positive function on 
an open neighborhood U of z e whose level sets are, by the Morse lemma, diffcomorphic to concentric 
spheres centered at the equilibrium z e . Additionally, both hypotheses (i) or (ii) imply that the constant 
e can be chosen small enough so that the time derivative 

L e (z) = {L e , H}(z) = - F(z e )) 2 , H}(z) + e^{F - F(z e ), H}(z) < 0, (2.3) 

for any z G P \ {z e }. This implies that if F t is the flow of Xh, the basis of open neighborhoods of 
z e given by the sets U\ := L7 1 ([0, A)), with A small enough, satisfies F t {U\) C F S (U\), provided that 
t > s. This proves the weak asymptotic stability of z e . 

If {(F — F(z e )) 2 , H}(z) < for every z £ P \ {z e } then e can be chosen so that the positive function 
L e is such that L e {z) < for any z G P \ {z e } (see H2.3fl ). Lyapunov's theorem proves the asymptotic 
stability of z e . ■ 

Remark 2.6 The main differences between this result and those already existing in the literature are: 

(i) It takes advantage of the possible existence of strict Lyapunov functions and hence is capable of 
obtaining the Lyapunov stability of an equilibrium as a corollary of an asymptotically stable 
behavior. This feature allows us to prove stability in some examples where no other available 
energy method is applicable. 

In order to illustrate this point consider the following example. The two dimensional Toda lattice 
admits a Poisson formulation |B100| by taking the bracket {x, y} = —x and the Hamiltonian 
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function H(x 1 y) — x 2 + y 2 . This system has an equilibrium point at z e = (0,6) for any del 
The equilibrium (0,0) is obviously Lyapunov stable since diJ(0, 0) = and d 2 iJ(0,0) > 0. The 
equilibria of the form z e = (0, b) with b > are weakly asymptotically stable. This can be 
proved using the previous theorem by taking the Hamiltonian as conserved quantity and the 
function F(x, y) :— x. The function F satisfies the hypothesis (i) in Theorem 12 . 51 since {F 2 , H} = 
—Ax 2 y < and {F, H} — —2xy — when {F 2 ,H} = 0, in an open neighborhood of z e = (0,6) 
with b > 0. If b < the equilibrium is unstable since the linearization has an eigenvalue with 
positive real part. We emphasize that the stability of the points in the case b > are uniquely 
due to their weak asymptotically stable behavior. 

(ii) Unlike the approach taken in the treatment of many standard examples (see for instance Ma Ra99] ^) 

this theorem shows that one does not need to take arbitrary functions of the conserved quantities 
in the expression 1)2.2(1 , Indeed, only linear combinations are needed. This is a consequence of the 
fact that the form whose definiteness needs to be studied is restricted to the space W. 

(iii) Since the constants {Ao, Ai, . . . , A&, fi} are allowed to be zero we have the freedom not to use a local 
conserved quantity in the definiteness condition (|2.2() but to still take advantage of its existence to 
shrink the space W. This is an improvement with respect to the results in |Or98| (see Corollary 
1-11 in |OrKa!)!)l,p 

In order to visualize this better consider the following example. Let (K 3 , {•,•}, H) be the Pois- 
son dynamical system whose Poisson bracket is given by the Poisson tensor that in Euclidean 
coordinates takes the form 

/ y\ 

B(x,y,z) = -a: 

\-y * o J 

and where H{x, y, z) — az with nela nonzero constant. The function C(x, y, z) — \ {x 2 + y 2 ) is 
a Casimir for this Poisson structure and every point of the form (0, 0, z ) is an equilibrium of Xh ■ 
Note that d(H — AC)(0, 0,z ) 7^ for any A e M. Nevertheless, we can still apply the previous 
theorem to conclude the Lyapunov stability of (0,0, zq) by taking the combination Ao-ff + AiC 
with Ao = and Ai = 1. With these choices, W = ker (dH(0, 0, z )) and d 2 C(0,0, z )\wxw 
is positive definite. The stability of these equilibria can also be handled using the topological 
methods in [Paal02j . 

Remark 2.7 The most efficient way to apply Theorem l2.5l in order to establish the stability of a given 
equilibrium consists of looking at the Hamiltonian system obtained by restriction of the original one 
to an arbitrarily small neighborhood of the equilibrium. The advantages of proceeding in this way are 
based on the fact that the restricted system has, in general, more conserved quantities than the original 
one. We illustrate this remark with the following specific example. 

Consider the manifold P := T 2 xl endowed with the Poisson structure given by the tensor that in 
coordinates (9, tp, x) is expressed as 

/0 1\ 

B(9,ip,x)=\ -a , ael\Q. 
\-l a / 

Let H € C°°(P) be the function defined by H(8,tp,x) := x 2 — cos 9. The associated Hamiltonian vector 
field Xh has an equilibrium at the point z e := (0, 0, 0) whose stability we show using Theorem l2.5l Even 
though the Poisson manifold P has no globally defined Casimir functions, any locally defined function 
of the form C — a9 + ip is a local Casimir. We can use this local Casimir to establish the Lyapunov 
stability of z e . Indeed, dH(z e ) — and d 2 H (z e )\wxw > 0, with W = kerdC(z e )- In section ET2l we 
will describe a mechanical system that is closely related to this example. 

Remark 2.8 In most applications, the conserved quantities in the statement of the theorem are local 
Casimir functions, components of momentum maps, and the Hamiltonian. A good way to find the 
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functions P is to look for purely negative eigenvalues of the linearization of Xh at the equilibrium z e 
that do not have a positive counterpart, as will be shown below. Notice that by Corollary 12.41 this is 
only possible when the equilibrium z e is lying on a singular symplectic leaf of the Poisson manifold. 
More explicitly, suppose that the linearization has such a negative eigenvalue —A with eigenvector v. 
Take local coordinates (yi, . . . , y n ) such that v = -Mj-. Since the function F v (yi, . . . , y n ) := y n satisfies 
{y^, H} — Xn[yn] — 2y n y n — —2\y^ + h.o.t., it is a good candidate to be used as the function F in the 
statement of the theorem. This procedure has been used in the first example in Remark 12.61 

2.3 Ideal reduction and ideal stability for Poisson systems 

We start this section by describing new Poisson structures on some submanifolds of a Poisson manifold 
that can be obtained by looking at the ideals of its Poisson algebra of smooth functions. We will refer 
to the construction that will be presented as ideal reduction for it is a particular case of the Poisson 
reduction procedures in |Mafia86l [OrRa98 , OrR a03j . 

This reduction technique is used later in this section to define a weaker notion of stability, called 
T-stability, and to establish a sufficient condition for it to hold. As the examples in the next section 
show, the use of I-stability is a very sensible way to deal with the physically relevant stability properties 
of equilibria in Hamiltonian systems subjected to semi-holonomic constraints. 

Let P be a smooth manifold and T C C°°(P) be a family of smooth functions. Denote by Vjf C P 
the vanishing subset of F, defined as the intersection of the zero level sets of all the elements of 
F. For a subset S C P define its vanishing ideal I{S) as the set of functions / € C°°(P) such that 
f(S) = {0}. Notice that I(S) is obviously an ideal of C°°(P) with respect to the standard multiplication 
of functions. Notice also that for every subset S C P and for every ideal J E C°°(P) we have S C Vi(s) 
and J Cl (Vj). These inclusions are in general strict. However, if S is a closed embedded submanifold 
of P then the first inclusion is actually an equality due to the smooth version of Urysohn's lemma. 
Moreover, in this particular case, the quotient algebra C°°(P)/X(S) can be identified with C°° (<!?), 
the algebra of smooth functions on S with respect to its own smooth manifold structure, via the map 
that assigns to any / e C°°{S) the element ir(F) 6 C°°(P)/T(S), where F £ C°°(P) is an arbitrary 
extension of / and tt : C°°(P) -» C°°{P)/1{S) is the projection. We will say that an ideal 1 C C°°(P) 
is regular if its vanishing set Vi C P is a closed and embedded submanifold of P. 

In the sequel we will focus our attention on finitely generated Poisson ideals. Let (P, {•,•}) be a 
Poisson manifold and F — {/i, . . . , /„} C C°°{P) be a finite family of elements in C°°(P). We will say 
that F generates a Poisson ideal if for any function / 6 C°°(P) and any i £ {1, . . . ,n} there exist 
functions {ha, . . . , hi n } C C°°(P) such that 



note that the condition above is equivalent to the statement that F(F) is an ideal in the Poisson algebra 
C°°(P), that is, it is an ideal relative to both the usual multiplication of functions as well as the the Lie 
bracket {•, •}. Note that if the vanishing subset Vjf of F is an embedded submanifold of P then Vjf is a 
quasi-Poisson submanifold of P. Indeed, for any / g C°°(P), fi 6 F, and z £ Vjf, there exist functions 



n 



{f,fi} = J2 h vh- 



3=1 



Denoting 




{hi,..., 



K] c C°°(P) such that 



11 



(dfi(z),X f (z)) = {fi,f}(z) =^2hi(z)Mz) = 0, 



which shows that B$(z)(T*P) C T z Vf, as required. Since the embedded submanifold Vf is quasi- 
Poisson, it has a Poisson bracket {-,-}v F given by {/, g}x> F (z) := {F,G}(z), where P, G £ C°°(P) are 
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arbitrary local extensions of /, g G C°°(Vf) around the point z eVf- We recall that the extensions to 
P of the Casimir functions of (Vf, {•, - }v F ) are called sub-Casimirs of (P, {-,-}). 

The construction that we just carried out can be locally reversed, that is, given an injectively 
immersed quasi-Poisson submanifold S of (P, {•, •}) any point z E S has an open neighborhood V z of z 
in S such that the vanishing ideal X(V Z ) is a Poisson ideal generated by a finite family of smooth functions 
on P with codimS' elements. Indeed, choose V z small enough so that it is an embedded submanifold of 
P and that, at the same time, is contained in the domain of a submanifold chart (U z , (f) of P. With this 
choice we can write U z ~ W\ x Wi and V z — Wi x {0}, where W\ and Wi are open neighborhoods of 
the origin in two finite dimensional vector spaces of dimensions dim S and codim S, respectively. If we 
denote the elements of Wi by (xi, . . . , a^codims) then any arbitrary extensions Pi, ... , P co dimS G C°°(P) 
of the coordinate functions fx — x%, . . . , / co dimS = ^codimS to the manifold P generate 1{V Z ) and form 
a Poisson ideal. Indeed, since V z is an embedded quasi-Poisson submanifold of P, we have for any 
F G C°°(P) and any z' G 7* 

{P,P}(z') = {/ 4 ,P|yJ^(^')=0 

since fi\v z = 0. 

Some of the ideas that we just introduced play a very important role in the algebraic approach to 
Poisson geometry. The reader interested in these kind of questions is encouraged to check with [Va96 
and references therein. 

Definition 2.9 Let (P, {•,•}, H) be a Poisson dynamical system and let T be a regular Poisson ideal, 
that is, the vanishing set Vz is a closed and embedded submanifold of P. Consider the reduced Poisson 
system (Vz, {•, -}vx, h) where h £ C°°(Vz) is defined by h := Hoi with i : Vj P the inclusion. Assume 
that z e G Vz is an equilibrium point for the Poisson dynamical system (P, {■,■}, H) and hence also for 
(Vr, {•) - }vx) h)- We say that z e G Vr C P is an X-stable equilibrium if any of the two following 
equivalent conditions hold: 

(i) z e is a stable equilibrium for the reduced Poisson dynamical system (Vz, {■, -}vx) h); 

(ii) for any open neighborhood U of z e in P , there is an open neighborhood V of z e in P such that if 

F t is the flow of Xu, then F t {z) £ U D Vz for any z G V R Vz and for any t > 0. 

The equilibrium z e is X-unstable if z e is an unstable equilibrium for the reduced Poisson dynamical 
system (Vz, {•, - }vx> ^ * s obvious that X -instability implies Lyapunov instability on the whole space. 

Theorem 2.10 Let (P, {•, •}, H) be a Poisson dynamical system with an equilibrium at the point z e G P 
and let U C P be an open neighborhood around z e . Assume that there exists a regular Poisson ideal 
X generated by the functions G\, . . . , G rn G C°°(P) with sub-Casimirs Pi, . . . , F r G C°°(P) such that 
z e G Vz- Suppose that the functions Co := H, Ci, . . . , C n G C°°(P) are conserved by the flow of Xr and 
that, additionally, there exist constants Ai, . . . , A n , ui, . . . , u r , v\, . . . , v m such that 

(i) Hi :— Y^7=o ^i^i + X^=i Mj'-^j' + SfcLi ^fcGfc /ias a critical point at z e , and 

(ii) i/ie Hessian of H\ at z e is positive definite when restricted to the subspace W defined by 

n r m 

W=f] kev(dCi{z e )) f] ker(dP J (z e )) f| ker(dG fc (z e )). 

i=0 j=l k=l 

Then z e is an X-stable equilibrium. 

Proof The hypotheses in the statement of the theorem imply that the equilibrium z e of the reduced 
system (Vz, {•, -}vn H o i), with i : Vz P the inclusion, satisfies the hypotheses of Theorem 12 . 51 and 
hence is Lyapunov stable on Vz, which implies that z e is I-stable. I. 
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3 Examples. 

3.1 A light Chaplygin sleigh on a cylinder. 

The following example was formulated in Mar95 in the context of nonholonomically constrained sys- 
tems. In that work the author found an equilibrium that exhibits asymptotically stable behavior. Wc 
will study the stability of all the equilibria of this system as well as of its relative equilibria with respect 
to a circle symmetry of the system that will be introduced later on. We will apply the Lyapunov sta- 
bility methods presented in the previous sections. This example is based on a real mechanical system 
that illustrates the theory particularly well since it exhibits equilibria that are not critical points of 
the Hamiltonian or of any other conserved function and, nevertheless, Theorem 12.51 still allows us to 
establish the Lyapunov stability of some dynamical elements and, in some cases, asymptotic stability. 
There is also an equilibrium to which none of the stability methods in the paper apply but that, after 
ideal reduction, is shown to be unstable and hence unstable in the whole space. 

We will start the presentation by explicitly carrying out in this particular example Marie's reduction 
procedure for nonholonomically constrained systems. The reader is encouraged to check with the original 
references |Mar95l Mar98 in order to find various technical details that we will omit here. 

3.1.1 Description of the system. 

The configuration space is given by the points (x,9) on a cylinder Q :— R x S . The Lagrangian 
of the system is just the kinetic energy L = \{x 2 + 2 ) £ C°°(TQ). The system is constrained to 
move subject to the semiholonomic constraint x + x9 = 0. The term "semiholonomic" means that 
the distribution that describes the constraint is integrable with integral leaves that are not necessarily 
embedded submanifolds. 

This system approximates a simple mechanical system in a certain regime that can be physically 
realized in the following way. Take a Chaplygin sleigh moving in the interior of a cylinder (we are 
assuming that all the physical constants of the system are equal to 1). The configuration space of 
this system consists of the points (x,8,tp) 6 Q' :— M x S 1 x S 1 , where the coordinates (x,8) on the 
cylinder indicate the position of the Chaplygin sleigh. The dynamics of this system is determined by 
the Lagrangian V on TQ 1 given by L' = \(x 2 + B + I^fi 2 ), where I v is the moment of inertia of 
the sleigh, together with the nonholonomic constraint x cos tp — 6 sin if = 0. Assume now that we add 
a new holonomic constraint tan<^ = x. Notice that even if the first constraint was not integrable, 
the superposition of the two constraints is integrable. In this case the dynamics can be described by 
restricting the system to a new configuration space Q C Q' which is actually an integral manifold of the 
distribution that describes the holonomic constraint. Moreover, it is easy to see that we can restrict the 
system to the integral manifold of any subset of integrable constraints, obtaining a new holonomically 
constrained system. In this case, we restrict the system described by the Lagrangian L' on TQ' to 
the integral submanifold Q C Q' by using the holonomic constraint tan ip = x. Assuming I v -C 1 
and restricting our study to points such that i C 1, the example that we will be presenting is a good 
approximation of this mechanical system. Marie |Mar95| considers the same mechanical realization of 
these equations but he sets I v = from the beginning of his exposition. 

3.1.2 Reduction of the system 

We now apply a reduction procedure due to Marie |Mar95l [Mar98 to eliminate the semiholonomic 
constraint x + x6 = 0. This reduction procedure consists of eliminating the Lagrange multipliers of a (in 
general nonholonomically) constrained system by finding a submanifold (the constraint submanifold) 
endowed with an almost Poisson structure and a Hamiltonian on it in such a way that the dynamics of 
this almost Poisson dynamical system coincides with the dynamics of the original constrained system. 
There are several equivalent constructions (see vdSMa94 Cual95, Mar95 Blal96, SniaOl , and references 
therein) to handle these constraints. It was shown in |vdSMa94] that this almost Poisson structure is 
actually Poisson if and only if the constraints are semiholonomic. 
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Let Q = M x S 1 be the configuration space and L(x, 9, x, 9) = \(x 2 +9 ) the Lagrangian of the system 
subjected to the constraint x + x9 = 0. Since the Lagrangian L is hyperregular, the Legendre transform 
FL : TQ — > T*Q is an isomorphism that we use to associate a Hamiltonian function H G C°°(T*Q) to 
the system. The image by FL of the constraint submanifold in TQ gives the constraint submanifold P 
on T*<3 which consists of the points P = {(x, 9,p x ,p ) G T*Q \ p x + xp = 0}. Let V C T{T*Q) be the 
so called constraint distribution defined by T>(z) := T Z P for every z in P. D'Alembert's principle 
provides a prescription to modify the original unconstrained Hamiltonian flow in order to construct 
a new vector field whose integral curves lie in P. Indeed, let Xh\ p be the restriction of the original 
Hamiltonian flow to the points in P and let Xu be the modified vector field whose integral curves describe 
the dynamics of the nonholonomically constrained system. The works by Marie quoted above ensure 
that, under certain regularity conditions satisfied in this example, the difference Xw = Xh\ p — Xjj of 
these two vector fields, is a section of a subbundle W of Tp(T*Q) that satisfies T P (T*Q) — W © V and 
that is uniquely determined by D'Alembert's principle. In such a situation, every Hamiltonian vector 
field can be decomposed in a unique way as Xh \ p = Xd + Xw and Xd describes the dynamics of the 
constrained system. Marie also shows that there exists an almost Poisson structure on P with almost 
Poisson tensor B : T*P x T*P — > R, for which X D = B*dH\ P , where B B : T*P -> TP is the canonical 
vector bundle isomorphism associated to B. 

In our example, V(x,9,p x ,p e ) = span{(l, 0, -pe, 0), (0, 1, 0, 0), (0, 0, -x, 1)} and W {x, 9,p x ,p e ) = 
span{(0,0, l,a:)}. An explicit expression for the almost Poisson structure (see |OrP104| ) can be given 
by using the natural projection map onto the T> factor. After some computations this almost Poisson 
tensor takes the form: 





B(x,0,p e )=\ 

x —1 



1+x 2 





1+x 2 

where the three-tuples (x, 9,pg) are used to coordinatizc the points (x, 9, —xpg,pg) € P and the restricted 
Hamiltonian is given by H\p(x, 9,p$) = |(1 + x 2 )pg. Notice that this tensor is Poisson since the 
constraint is integrable. The equations of motion are 

x = -xp e , 9 = po , f, 



1 + x 2 



3.1.3 Equilibria, relative equilibria, and their stability 

Notice that every point of the form z = (x, 9, 0) is an equilibrium of the system. If we first compute the 
linearization of the dynamical system at those equilibria we obtain the family of matrices 




which have three zero eigenvalues and are not diagonalizable. This implies that the system is linearly 
unstable at those equilibria (which does not imply either Lyapunov stability or instability). 

To apply Theorem 12.51 we first need to find conserved quantities for the Hamiltonian flow. In this 
case we can use the Hamiltonian and the local Casimir function given by C(x,9,pg) — xe e . Let L be 
the function defined by L := XqH + \\C . If we set Ao = 1 and Ai = we have that dL(z) — 0. The 
subspace W — kerdH(z) n dC(z) is given by W = span{(x, —1,0), (0,0, 1)} and the restricted Hessian 





d2L W\wxW-y (l + x 2 ) 

is not positive definite since it has a zero eigenvalue. The stability of the equilibrium z = (0, 9, 0) can be 
analyzed by using the fact that the submanifold S consisting of the points of the form (0, 9,pg) is such 
that its vanishing ideal T{S) is a Poisson ideal and hence S is Poisson reducible. Indeed, if (9 7 pg) are 
coordinates on S, the reduced bracket {•, -}g takes the form {9,pg}s = 1 and the reduced Hamiltonian 
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is h(9,pg) — hp$- This reduced system describes a free one dimensional particle. The equilibrium 
z = (0, 8, 0) of the original system drops to an equilibrium at the point (8, 0) which is clearly unstable. 
In particular, this implies the instability of the original equilibrium (0, 8, 0). 

We now study the stability of the relative equilibria with respect to the circle symmetry of the 
system given by the action tp ■ {x,9,pg) = {x, 9 + ^,pg). This action is canonical and the system can 
be Poisson reduced. The reduced manifold is R 2 . If we denote by (x.pg) the elements of the reduced 
space, the reduced Poisson bracket is determined by the relation {x, pe} — —x/{l + x 2 ) and the reduced 
Hamiltonian is h(x,pe) = i(l + x 2 )p 2 e . Hamilton's equations for h are x — —xpg, pg = x 2 pg/(l + x 2 ). 
Thus the equilibria are given by the family of points satisfying xpg = 0. The linearization of the 
Hamiltonian vector field at these equilibria is given by the matrix 



which has a positive eigenvalue if pg < 0, in which case the system is Lyapunov unstable at the 
points (0,pe). This obviously implies that the unreduced system exhibits nonlincarly unstable relative 
equilibria. 

If pg > the linearization does not imply neither stability nor instability. However, note that in 
this case, the linearization has a negative eigenvalue with eigenvector v — (1,0) that will be useful 
when searching for a Lyapunov function (see Remark 12.8(1 . In order to study the nonlinear stability of 
these relative equilibria, we notice that the only available conserved quantity is the reduced Hamiltonian 
whose derivative dh{x,pg) — (%Pg, (1 + x 2 )pg) = (0, 0) if and only if pg = 0. In that case 



and hence we cannot conclude either stability or instability. However, in this particular case instability 
can be concluded just by looking at the phase portrait for the vector field. For points of the form (0,pg) 
the derivative of the Hamiltonian does not vanish and hence the only way to apply Theorem l2 . 51 consists 
of finding a function F satisfying at least one of the hypotheses (i) or (ii); F(x,pg) = x 2 /2 is one such 
function since {x 2 , h} = —2x 2 pg, {x 4 , h} = —4x 4 pg, and pg is assumed to be positive. Consequently, the 
hypothesis (i) is obviously satisfied. With this choice, the subspace W = ker dh(0,pg) = span{(l,0)} 
and d 2 F(0,pg)\ Wxw — 1 > 0. Consequently, the equilibria of the form (0,pg) with pg > are Lyapunov 
stable and even though they are not asymptotically stable, there exists an open neighborhood V of (0,pg) 
such that F t (V) C F S (V), whenever t > s, that is, they are weakly asymptotically stable. 

Finally, it is easy to conclude that the equilibria on the form {x, 0) are unstable just by looking at 
the phase portrait of the system. 

3.2 Two coupled spinning wheels 

Consider two vertical weightless wheels with radii R and r satisfying R > r and R/r G R\Q. We attach 
to the edges of each of these wheels two point masses M and m. This simple system has as configuration 
space Q the torus T 2 that we coordinatize with the angles (8, ip). The Lagrangian of this system in these 
coordinates is L = ^(MR 2 9 +mr 2 Lp 2 ) + AIR cos 8 + mr cos ip. Assume now that we couple the rotations 
of the two wheels with a belt. This mechanism imposes on the systems a semiholonomic constraint that 
can be expressed as R9 — rip = 0. In order to give a description of the constrained system we first 
express the original system in the Hamiltonian setting by using the Legendre transform. The phase 
space P is in this case the cotangent bundle T*T 2 ~ T 2 x R 2 with coordinates (8,ip,pg,p v ), endowed 
with the canonical symplectic form. The Hamiltonian function is 



The constraint submanifold is given by the points(#, ip,pg,p v ) that satisfy p v = mrpg/MR, which can 
be identified with T 2 x R with coordinates (9, ip,p). 
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We now apply the reduction procedure in |BaSn93] in order to find a bracket on the constraint 
submanifold that is actually Poisson since the constraint is semiholonomic. This bracket is given by the 
constant Poisson tensor: 

/ r\ 
B(9,<p,p) = R). 

\-r -R 0/ 

The reduced Hamiltonian function is 

P 2 

h(0, <p,p) = — — MRcosd — mr cos <p, 
2k 

where k is a real positive constant depending on the parameters of the problem given by the expression 

m + M (to - M) 2 m 2 M 2 

= 4MmR 2 r 2 ~ 4m 2 M 2 R 2 r 2 (m + M) ' 

This Poisson system has a local Casimir given by the locally defined function C(9, ip,p) = R0 — rip. The 
equations of motion of the system are given by 

P.p. 
6 = r—, (p — R—, p = —rRM sin 9 — mrR sin w. 
k k 

The equilibria of the system are the points of the set S = {(6, <p, 0) | Msin6> + to sin <p = 0} that 
can be described as a one-parameter family given by the curve (p = — sin -1 ( M T s " e ) , G [— 9 C , 6 C ], 
where 9 C is given by 9 C = sin (■?£■). In order to study the nonlinear stability of such equilibria we 
compute Xodh(z) + XidC(z) = 0, with z = (6, ip, 0) e S. This equation can be solved by taking 
Ai = -X M sin9. In this case W — kerdC(z) nkerdft-(z) = span{(r, R, 0), (0, 0, 1)}. Finally, it is easy 
to see that d 2 (Xoh + XiC)(z)\wxw > if and only if Mr cos 9 + mR cos ip > 0. In particular, the point 
z — (0, 0, 0) is always nonlinearly stable, as expected, and the point z = (0, 7T, 0) is stable if ^ > y . 

4 Nonlinear stability via topological methods 

In |Paal02| topology based tools have been developed that provide sufficient conditions for the Lyapunov 
stability of Poisson equilibria. One of the main achievements in Pa al02) is the discovery of a space 
related to the topology of the symplectic foliation of the Poisson manifold (see l|4.3|l below) on which 
the extremality of the Hamiltonian suffices to conclude stability. In this section we will study under which 
circumstances the topological criteria in |Paal02| can be expressed in terms of local continuous Casimir 
functions and hence there is an equivalence with the energy-Casimir method. To be more explicit, we 
will seek the correspondence between the topological approach of Paal02j and a generalization of the 
energy-Casimir method that requires only continuity of the functions involved and that is based on the 
following general lemma. 

Lemma 4.1 Let X G X(P) be a smooth vector field on the finite dimensional manifold P and z e G P 
an equilibrium point. If there exists locally defined continuous conserved quantities Co, . . . , CV- G C (U) 
of the flow F t such that f] i=0 C~ 1 (Ci(z e )) — {z e } then the equilibrium z e is Lyapunov stable. 
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Proof Consider the function L(z) = {C (z) — Co(z e )) 2 + ... + (Cfc(z) - C'k(z e )) 2 . The hypothesis 
HiLo C~ 1 (Ci(z e )) — {z e } ensures that L is a positive function that takes the zero value only at the point 
z e . In particular, the sets of the form i _1 ([0, e)), e > 0, form a fundamental system of neighborhoods 
in the manifold topology of P at the point z e . Consequently, for any open neighborhood U of z e there 
exists an e > such that L _1 ([0, e]) C U. Since the level set L~ 1 (e) is invariant by the flow F t of X, 
the Lyapunov stability of z e follows. I 

Any continuous function C G C (U), with U an open subset of P, such that C is constant on 
the symplectic leaves of (U, {•, -}\u) is called a local continuous Casimir of (P, {•,•})■ The choice 
of terminology is justified by the fact that if such a function C happens to be differentiable then it 
is an actual Casimir of (U, {•, -}\u)- It is worth noticing that the local continuous Casimirs are the 
(continuous) first integrals of the foliation of (U, {•, -}\u) by its symplectic leaves. 

Corollary 4.2 (Continuous energy-Casimir method) Let (P, {•,•}, H) be a Poisson dynamical sys- 
tem and z e 6 P an equilibrium point of the Hamiltonian vector field Xh ■ Let S Ze C P be the common 
level set of local continuous Casimir functions around z e . If 

H-\H(z e ))nS Zc ={z e } (4.1) 

then the equilibrium z e is Lyapunov stable. This statement remains true if H is replaced by any contin- 
uous conserved quantity of the flow of Xh . 

Our goal is to establish sufficient conditions under which this corollary coincides with the topological 
stability criterion in Paal02 that we now recall. We start by introducing the necessary notation. 
Let (X,t) be a topological space and x G X an arbitrary point. We define the set T 2 (x) clas 

T 2 (x) := {y G X | U x fl U y ^ for any two open neighborhoods U x , U y of x and y}. (4.2) 

Let A C X be an arbitrary subset. We define 

T 2 (A) := |J T 2 (x). 

Notice that if y G T 2 (x) then x £ T 2 (y). Also, a topological space (X,t) is Hausdorff if and only 
if T 2 (x) — x for every x S X. Hence the T 2 sets measure how far a topological space is from being 
Hausdorff. 

Suppose now that P is a smooth Hausdorff and paracompact finite dimensional manifold and D is 
a smooth and integrable generalized distribution on P. Let 7Td : P — » P/D be the projection onto the 
leaf space of the distribution D. The map n is continuous and open when P/D is endowed with the 
quotient topology. Define 

T 2 (x)=TTp 1 (T 2 (ir D (x))), xeP (4.3) 

and, more generally, 

T 2 (x)=n D l(T 2 (7r Dlu (x))), x E P, (4.4) 

where U is an open neighborhood of x 6 P and i^d\u '■ U — > U/D\u is the projection onto the leaf space 
of the restriction D\u of D to U. 

We now focus on the particular case when P is a Poisson manifold with bracket {•,•}■ Let £ be 
the corresponding characteristic distribution and ir : P — > P /{■,■} the projection onto the space of 
symplectic leaves P/{-, ■} := P/£. 

Theorem 4.3 (Topological energy-Casimir method [Paal02 ) Let (P, {•,•}, H) be a Poisson dy- 
namical system and z e G P an equilibrium point for the Hamiltonian vector field Xh ■ If there is an 
open neighborhood U C P of z e such that 

H- l {H{z e ))^ 7 f J 2 {z e ) = {z e } (4.5) 

then the equilibrium z e is Lyapunov stable. This statement remains true if H is replaced by any contin- 
uous conserved quantity of the flow of Xh that takes values in a Hausdorff space. 
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In view of expressions (|4.1() and l|4.5|l we would like to know under what circumstances the set T 2 (z e ) 
can be obtained by looking at the level sets of local continuous Casimir functions thereby rendering the 
statements of Corollary 14.21 and Theorem ED equivalent. 

The first point that we have to emphasize is that this is, in general, not possible. The following 
example, that we owe to James Montaldi, shows that, in general, we cannot find enough local Casimir 
functions to be able to write the set T 2 (z e ) as the common level set of local continuous Casimir functions, 
no matter how much we shrink the neighborhood U. Let K 3 and f(x,y,z) = x 2 + y 2 — z 2 . Consider 
the Poisson structure {•, •} determined by {x, y} — f 2 , {y, z} = 2yzf, and {x, z} — —2xzf. In order to 
describe the symplcctic leaves of (R 3 , {•, ■}) (see Figure}? jj notice first that the function / is a factor and 
hence the Poisson tensor vanishes on the cone / = 0. Consider now all the spheres through the origin and 
tangent to the OXY plane (and hence centered on the OZ-axis) and cut them with the cone / = 0. Each 
of these spheres contains the following symplectic leaves: the sphere intersected with the points (x, y, z) 
such that f(x, y, z) > (two dimensional leaf), the sphere intersected with the points (x, y, z) such that 
f{x, y, z) < (two dimensional leaf), and the points such that f(x, y, z) = (zero dimensional leaves). 
It is clear from this description that there are no non-constant continuous local Casimir functions near 
the origin. Nevertheless, for any neighborhood U of the origin T 2 (0,0,0) = {(x,y,z) | f(x,y,z) > 0}, 
that is, the closed exterior of the cone, which in this case is strictly included in C^ l {C v {Q, 0, 0)) = U. 



Figure 4.1: Symplectic leaves of Montaldi's example of a Poisson manifold that does not have local 
Casimirs around the origin. The shadowed area represents the set T 2 (0,0,0). The picture is a section 
of the three dimensional figure through the OYZ plane. 

Even though the previous example shows that the set T 2 (z e ) does not coincide in general with the 
common level set of local continuous Casimir functions one can easily prove that at least one inclusion 
holds true. The natural context to present most of the results in this section is that of generalized 
foliations of smooth manifolds. Consequently, we will prove our statements in that category and we 
will obtain the Poisson case as a corollary by applying the theorems to the generalized foliation of the 
Poisson manifold by its symplectic leaves. 

Lemma 4.4 Let P be a smooth finite dimensional manifold and D a smooth integrable generalized 
distribution on P. Let ttjj : P — » P/D be the projection onto the leaf space of the distribution D and T 2 
the symbol defined in (|4.3(l . Let Ci G C°(P), i G I, be a set of continuous functions that are constant 
on the integral leaves of D (that is, first integrals of D). Then for any z G P 



Proof. Let C : P -> R 1 be the function defined by C{z) := (C t (z)) teI . If we endow M. 1 with the 
product topology (not the box topology!) then the continuity of the first integrals Ci, i G /, implies 
that C is continuous. The projection tt : P — > P/D is an open map when P/D is endowed with the 





(4.6) 
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quotient topology. Given that C is constant on the integral leaves of D it drops to a map c : P/D — > M. 1 
that closes the diagram 

C 

P - R 7 




P/D 

The continuity of C and the openness and surjectivity of ttd imply that c is also continuous. In order 
to prove H4.6f> it suffices to show that if m £ T2 (z) then C(m) — C(z). By contradiction, suppose 
that C(m) 5* C(z). Since R 7 is a Hausdorff topological space there are open neighborhoods Vc( m ) and 
^c(z) °f TO and z, respectively, such that Vc( m ) H Vc^z) =0. As c is continuous the sets c (Vc(m)) 
and c~ 1 (Vc( z )) are open neighborhoods of 7T£>(m) and 7Td(z)) respectively. Also, since by hypothesis 
m G T 2 (z), we have that c _1 (Vc( m )) n c~ 1 {Vc( yZ )) $j necessarily. However, by construction we also 
have that c~ 1 (V C ( m )) n c _1 (V C ( z )) = c _1 (Vc( m ) ("1 Vcr z )) = c _1 (0) = 0, which is a contradiction. ■ 

The rest of this section is dedicated to the description of two situations where the inclusion (|4.6I) 
is an equality and hence local continuous Casimir functions characterize the T2-sets. We start with a 
couple of preliminary general results. 

Definition 4.5 Let (A, r) be a topological space. We say that (A, r) is T 2 -idempotent when T^^Tiix)) = 
T^ix), for any x G X . 

Lemma 4.6 Let (X,t) be a T^-idempotent topological space. Then 

(i) The relation 1Zt 2 on X defined by xR-T 2 y if and only if y G T^(x) is an equivalence relation on X. 

(ii) The following statements are equivalent: 

1. y£T 2 (x). 

2. T 2 (x)^T 2 (y). 

3. T 2 (x)nT 2 (y) =0. 

4- There exist open neighborhoods U x , U y of x andy, respectively, such thatT 2 (U x )nT 2 (U y ) = 0. 

(iii) If the projection ttt 2 '■ A — ► X/TZt 2 onto the space of equivalence classes endowed with the quotient 
topology is an open map then X/1Zt 2 is a Hausdorff topological space. 

Proof (i) The definition of the T 2 set implies that x!Zt 2 x for any x G A and that xTZt 2 V if and only 
if yR-T 2 x. In order to prove transitivity of 7&t 2 let x,y,z(^Xbe such that xTZt 2 V and yTZx 2 z. By the 
very definition of the T 2 set, it is clear that for any two subsets A, B C A such that icBwe have 
that T 2 (A) C T 2 (B), in particular, the condition x G T 2 (y) implies that T 2 (x) C T 2 (T 2 (y)) = T 2 (y). By 
reflexivity we have that T 2 (y) C T 2 {x) and hence T 2 {x) = T 2 {y) which implies that T 2 {x) = T 2 (y) = 
T 2 {z) and hence x!Zt 2 z. 

(ii) If T 2 (x) = T 2 (y) then y G T 2 (y) = T 2 (x). This proves the implication 1=>2. The implication 2=>1 
was already proved in the first part of the lemma. In order to prove 2=^3 suppose that there exists a 
point z G T 2 (x) fl T 2 (y). Then using the T 2 idempotency as we did in the proof of the first part of the 
lemma we obtain that T 2 (x) = T 2 (z) = T 2 (y), which contradicts the hypothesis. To show 3=>4, assume 
that T 2 (x) n T 2 (y) = 0. Then, in particular, y T 2 (x) and hence there exist open neighborhoods U x 
and U y of x and y, respectively, such that U x fl U y = 0. Since U x and U y are open neighborhoods of 
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each of their points, it follows that for every a x G U x and a y G U y the element a x £ T 2 (a y ). Using the 
implication 1=>3 that we have already proved, this shows that T2(a x ) (~l T2(a y ) = and hence 



(J TMjf) ( U T 2K)j = (J (T 2 (a x )f}T 2 (a y ) 



T 2 {U x )nT 2 {U y ) = 
Finally, the implication 4=^2 is straightforward. 

(iii) Notice first that for every subset A C X, we have that 7Ty 1 (7rT , 2 (A)) = T2(A). Let p,a G X/7Zt 2 be 
two points such that p ^ o and let x and y be two points in X such that ttt 2 ( x ) = /° an d 7r(y) = o". Since 
T 2 (x) 7^ T2 (y) there exist, by part (ii), two open neighborhoods V x and V y of a; and y, respectively, such 
that = T 2 {V X ) n T 2 (V y ) = n£fa- a (V x )) D t^ttt^)) - ^Vr, (V K ) n ttt 2 (^)). Applying tt T2 to 
both sides of this equality we obtain that ttt 2 (Vx) H 7Tt 2 (Vy) = 0. Since ttt 2 (V x ) and 7Tt 2 (Vy) are, by the 
openness of ttt 2 , open neighborhoods of the points p and a, respectively, the claim follows. ■ 

Suppose now that P is a smooth Hausdorff and paracompact finite dimensional manifold and D is 
a smooth and integrable generalized distribution on P. Let 7r/j : P — > _P/L> be the projection onto the 
leaf space of the distribution D and T2 the symbol defined in 1|4.3JI . Notice that since 7td is surjective, 
we have 

ttd(T 2 (x)) = T 2 (n D (x)) , for any x G P. (4.7) 

We will say that the pair (P,D) is T 2-idempotent when T 2 (T2(x)) = 7 1 2 (a;), for any x £ P. Notice 
that since the sets T 2 (x) are D-saturated (they are unions of leaves of D) , we can conclude using (14.71) 
that P is T2-idempotent if and only if P/D is T2-idempotent. With this remark in mind the previous 
lemma can be easily adapted to the symbol T 2 . 

Lemma 4.7 Let P be a smooth Hausdorff paracompact finite dimensional manifold and D a smooth 
integrable generalized distribution on P. Let 7T/j : P — > P/D be the projection onto the leaf space of the 
distribution D and T2 the symbol defined in (|4.3(l . Suppose that (P, D) is T 2-idempotent. Then: 

(i) The relation Thp on P defined by xTZj^^y if and only if y £ T 2 {x) is an equivalence relation. 

(ii) The following properties are equivalent: 

1. y_i T 2 {x)_ 

2. T 2 (x) ^T 2 (y). 

3. T2(x)nT 2 (y) = 0. 

4- There exist open neighborhoods V x ,V y of x andy, respectively, such that T2{V x )nT2(V y ) = 0. 

(iii) If the projection Ttj; : P — > P/TZ?p^ is an open map then the quotient space P/TZ?p 2 is a Hausdorff 
topological space. 

Proof (i) Only transitivity needs to be proved. Let x,y,z G P be such that x12*p y and yTZ^^z. By 
definition, 7T/j(x)7?/r 2 7r.D(y) and i:D(y)'T^T 2 '^D(z). Since the T2-idempotency of (P,D) is equivalent to 
the T2-idempotency of P/D, Lemma FOl guarantees that , kd(x)TZt 2 '^d(z) and hence 7t_d(x) £ T 2 (ttd(z)). 
Consequently, x G tt/j 1 (T 2 (ttd(z))) — T 2 (z) and thus zlZ^x. 

In order to prove parts (ii) and (iii) it suffices to mimic the corresponding implications in Lcmma l4.6l 
but, this time, keeping in mind that the projection 7r^ 2 : P — > P/TZip^, = 7Tt 2 ° 7Td, is just the 
composition of two projection maps and that ttd is an open map. I 

Theorem 4.8 Let P be a smooth Hausdorff paracompact finite dimensional manifold and D a smooth 
integrable generalized distribution on P. Let 7T/j : P — > P/D be the projection onto the leaf space of the 
distribution D and T 2 the symbol defined in H4.3|) . Suppose that (P,D) is T 2-idempotent and that ttt 2 
(and hence ^x 2 ) * s °P en - Then the continuous first integrals of D separate the T 2 sets. In this situation, 
for any z G P , there exist continuous first integrals {C^jig/ C C°(P) of D such that 

T2(z) = f)Cr\C. l (z)). (4.8) 

iel 
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Proof Since P is by hypothesis paracompact, so are the quotient spaces P/D and P/IZ^^. The hy- 
pothesis on the T2-idempotency of (P, D) implies, by Lemma 14*771 that the quotient space P/Tl^ is also 
Hausdorff. Since a Hausdorff paracompact space is normal, Urysohn's Lemma guarantees the existence 
of continuous functions / on P/T?.^ that separate its points. The pull back / o Wj: G C°(P) is a first 
integral of D. The family of functions of the form foiTTp^ where / : P/lZj;^ — > R is a continuous function 
that separates two arbitrary points, is the family of continuous first integrals of D in the statement of 
the theorem. 

In order to prove the identity (|4.8|l it suffices to reproduce the proof of Lemma T4. 41 taking this time 
the function C : P — > R 7 whose components are the continuous first integrals of D that separate the T 2 
sets and whose existence we just proved. I 

Remark 4.9 The two hypotheses in the statement of this result, that is the T2-idempotency and the 
openness of the projection 7Ty 2 are independent. Indeed, consider the foliation of the Euclidean plane 
R 2 by the integral curves of the vector field <fi(x)d/dx, where is a smooth function satisfying <j)(x) = 0, 
for x < 0, and 4>{x) > 0, for x > 0. In this situation T 2 {x,y) = {(x,y)}, when x < 0, and T 2 (x,y) — 
{{ x ilj) G R 2 I 1 < 0}, if x > 0. In this situation, we obviously have T2-idempotency. However, the 
projection 7^ : R 2 — > R 2 /7^ 2 is not open. Indeed, the saturation n^^^^U)) = {(x,y) G R 2 | x > 0} 

of the open set U = {(x, y) G R 2 | x > 0} is closed, which is not compatible with 7Ty 2 being open. 

The following result provides another sufficient condition for the conclusion of Theorem l4.8l to hold. 

Theorem 4.10 Let D be a generalized smooth integrable distribution defined on the second countable 
finite dimensional manifold P. Suppose that there exist continuous first integrals Ci G C°(P), i G / , of 
the foliation induced by D that separate its regular leaves. Additionally, assume that the map C : P — > M. 1 
defined by C(z) := (Ci(z))i e j , z G P , is open onto its image when M. 1 is endowed with the product 
topology. Then for any z G P 

T 2 (z) = f]Cr 1 (C i (z)). 
iei 

Proof Notice first that the inclusion 

T 2 (z)c[]Cr 1 (C i (z)). 
iei 

is a particular case of H4.6|l . 

In order to prove the converse inclusion let ttd '■ P — * P/D be the projection onto the leaf space 
and c : P/D — > R fc the continuous mapping uniquely determined by the relation c o n D = C. Let 
n G C\ ieI C~ 1 (Ci(z)) , that is, C(n) = C(z) and assume that n £ T 2 (z). This implies the existence of two 
open neighborhoods V„ D r n ^ and V nn r z ^ of 7Td(«) and ttd(z), respectively, such that V nD t n ) n V 1TD ^ = 0. 
We will assume for the time being that the leaf 7r£)(n) is regular and will prove that the assumption 
n (/ T 2 {z) leads to a contradiction. We will prove later on that the situation in which ttd(ji) is a singular 
leaf can be reduced to this case. 

If 7ro(n) is regular, the set V^ g / n ) of regular leaves in V VD (n) is an open dense neighborhood of nr)(n) 
in V^ D (n)- The openness hypothesis on the map C implies that the set 

u c{z) :=c(q n) )nc(i; DW ) 

is an open neighborhood of C(z). Moreover, the continuity of c implies that the sets 
^:=c- 1 ((7 c(2) )n^„ ) and B := ^{U^)) n K D W 

are open neighborhoods of 7Td(ti) and ttd(z), respectively. Let ttd(z') be a regular leaf in B. The 
construction of B implies that there exists a regular leaf ttd(s) G V^^ 3 r n \ C Kr D (n) sucn that c(ttd(z')) = 
c(ttd{s)). The separation hypothesis on the map C implies that ttd(s) = no(z') G V^ D ( n ) C\V^ D ^ which 
is a contradiction. 
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In order to conclude the proof we need to show that the case in which 7ro(n) is singular can be 
reduced to the situation that we just treated. Indeed, take Uc{z) '■— C (K- D (n)) H c(V 7TD ^ z - ) ). By the 
openness of C, Uc( z ) is an open neighborhood of C(z). Additionally, the continuity of c implies that the 
sets A := cr 1 {Uc( z )) H K- D („) and B := c^ 1 (Uc( z )) H V-k d ( z ) ar e open disjoint neighborhoods of 7ro(n) 
and tto{z), respectively. Let 7Td(z') be a regular leaf in A. The construction of A implies the existence of 
a leaf itd(s) G V WD r z \ such that C(itd(z')) = C(ttd(s)). If we follow the preceding argument, replacing 
ttd(z') by 7T£i(n), 7T£>(s) by 7r/j(z), A by l / 7rD („), and with V^ D ( Z \ playing the same role we also obtain 
a contradiction with the hypothesis V^ D i n ) fl K- D (^) = 0- B 

The reader may be wondering how the two sufficient conditions for (|4.8() to hold that we presented 
in the statements of the theorems 14 . 81 and 14.101 are related. Our next result answers this question. 

Proposition 4.11 Let P be a smooth second countable finite dimensional manifold and D a smooth 
integrable generalized distribution on P. Suppose that there exist continuous first integrals Ci G C°(P), 
i G I, of the foliation induced by D that separate its regular leaves such that the map C : P — > C (P) C K / 
defined by C{z) := (Ci(z))i e i , z 6 P, is open onto its image when M. 1 is endowed with the product 
topology. Then (P,D) is 1 1 2 -idempotent and n-^ '■ P — ► P/TZ^ is an open map. 

Proof. In the hypotheses of the statement, Theorem 14.101 implies that T 2 (z) — C _1 (C(z)), for any 
zeP. In particular 

T 2 (T 2 (z)) = T 2 {C- l {C{z))) = |J T 2 (y) = C-\C{z)) - T 2 (z), 

yec-Hdz)) 

which guarantees that (P, D) is T2-idcmpotent and hence allows us to define an equivalence relation 'R^ 2 
on P. We will now show that the associated projection to the quotient 7^ : P — > P '/7&f is open. Let 
(p : P/TZ T2 -> C{P) be the map defined by (p(w T2 (z)) := C(z), z G P. The equality T 2 {z) = C-^C^)), 
z G P, guarantees that y is a well defined bijection that makes the diagram 

C 

P -C(P) 




Wt 2 

commutative. The continuity and the openness of C imply respectively the continuity and the openness 
of ip, that is, tp is a homeomorphism. Since 7^ = ip^ 1 o C, the openness of 7^ follows. I 

Remark 4.12 The converse of the implication in the previous proposition is not true in general. A 
counterexample to this effect is an irrational foliation of the two-torus. In that particular case the T 2 set 
of any point is the entire torus and hence we have T2~idempotency with a projection 7r^ : P — > P/TZ-^ 
that is obviously open. Nevertheless, the only first integrals of this foliation are the constant functions 
that do not separate the leaves of the foliation, all of which happen to be regular in this case. 

We now collect the results in Theorems 14 . 81 and 14. 101 and in Proposition 14. 1 II and we apply them to 
the situation in which P is a Poisson manifold foliated by its symplectic leaves. The following result 
provides two sufficient conditions for the continuous and topological energy-Casimir methods to coincide. 

Theorem 4.13 Let (P, {•,•}) be a Poisson (paracompact, second countable, and Hausdorff) manifold. 
Let T 2 be the symbol associated to the symplectic foliation of P induced by the Poisson structure {•, •}. 

(i) Suppose that there exist continuous Casimir functions Ci G C°(P), i G I, that separate the regular 
symplectic leaves of P such that the map C : P — > C(P) C M. 1 defined by C(z) :— (Ci(z))i^i , z G P, 
is open onto its image when M. 1 is endowed with the product topology. Then P is T 2 -idempotent 
and 7T7p 2 : P — ► P/7?.y 2 is an open map. 
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(ii) If (P, {•,•}) is T2-idempotent and 7Ty 2 : P — > P/TZj; is an open map then there exist continuous 
Casimir functions {Ci}i^i C C°(P) 0/ (P, {•, •}) swc/i £/ia£ /or any z G P 

T 2 (^) = ncr 1 (c i (^)). 

Remark 4.14 As one could expect, the hypotheses of this theorem are not satisfied by Montaldi's 
example (see Figure l4~T)l . Indeed, in this particular case T 2 (T 2 (0,0,0)) = U ^ T 2 (0,0,0), for any 
open neighborhood U of the origin (0, 0, 0). 

Example 4.15 The following example is given in |Paal02| and provides a situation where their topolog- 
ical energy-Casimir method fTheorem l4.3[) works when establishing the stability of a Poisson equilibrium 
while the standard energy-Casimir tool fails. We will show that in this situation the generalized energy- 
Casimir method formulated in Corollary 14.21 also works and, moreover, both results can be applied 
interchangeably since the hypotheses of part (i) in Theorem 14 . 1 31 are satisfied. 
Let (R 3 , {■,■}, h) be the Poisson dynamical system given by 

{/, g} = VA ■ (V/ x Vg), A(x, y, z) = (a 2 x 2 - y 2 )y 

where a is a nonzero real constant and h{x, y, z) — x 2 — y 2 + z 2 . Notice that the function A is a Casimir 
of the bracket {■, •} and that the points on the form (x,Q, z) are equilibria of the Hamiltonian vector 
field Xh- We will focus on the stability of the origin (0,0,0) that happens to be a singular point of 
the symplectic foliation of R 3 . In order to verify that the hypotheses of Corollary 14.21 are satisfied 
notice that the map A can be rewritten as A(x,y, z) = {ax + y)(ax — y)y and hence its zero level set 
(the one containing the equilibrium (0,0,0)) can be written as the union of three irreducible algebraic 
varieties V\, V 2 , and V3 that are the zero level sets of the functions y, ax — y, and ax + y, respectively. 
Consequently, 

h~ l (o) n A-\o) = ir\o) p) (Vi u v 2 u v 3 ) = (/^(o) n Vi) u (/^(o) n v 2 ) u (/^(o) n v 3 ) (4.9) 

which is a single point whenever \a\ < 1 hence proving the Lyapunov stability of (0,0,0). This is so 
since each of the three intersections on the right hand side of expression (|4.9|l coincide with the point m. 
This statement can be proved by showing that the Hamiltonian restricted to the submanifolds V\ , V 2 
and V3 has a non degenerate critical point at (0, 0, 0). This is closely related to the smoothing of the 
T 2 set introduced in IPaal02j . 

Since the Casimir function A clearly separates the regular symplectic leaves of (R 3 , {•,•}) and it is 
an open map, by Theorem 14 . 1 31 we can conclude that 

T 2 (0,0,0) = A- 1 (A(0,0,0)) 

and hence energy-Casimir and p2-based sufficient stability conditions can be used interchangeably. 
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